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Stability Analysis of Discrete-time Lur’e Systems with Slope-restricted 

Odd Monotonic Nonlinearities 
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Abstract — Many nonlinear dynamical systems can be written 
as Lur’e systems, which are described by a linear time-invariant 
system interconnected with a diagonal static sector-bounded 
nonlinearity. Sufficient conditions are derived for the global 
asymptotic stability analysis of discrete-time Lur’e systems 
in which the nonlinearities have restricted slope and/or are 
odd, which is the usual case in real applications. A Lur’e- 
Postnikov-type Lyapunov function is proposed that is used to 
derive sufficient analysis conditions in terms of linear matrix 
inequalities (LMIs). The derived stability critera are provably 
less conservative than criteria published in the literature, 
with numerical examples indicating that conservatism can be 
reduced by orders of magnitude. 

I. Introduction 

Absolute stability theory considers the stability of a nom¬ 
inal linear time-invariant system interconnected with a static 
nonlinearity, that is, a real function of real variables with 
no internal states. The literature on absolute stability theory 
can be categorized in terms of the properties of the static 
nonlinearities that are considered. This paper derives new 
sufficient conditions for the stability analysis of discrete-time 
systems with the advance being LMI conditions that exploit 
the static, sector-bounded, and slope-restricted nature of the 
componentwise nonlinear operators. 

Lyapunov methods provide a simple but powerful way to 
analyze nonlinear dynamical systems and design stabilizing 
controllers [1], [3], [14], [16], [35]. Many well-known re¬ 
sults on absolute stability were developed for a benchmark 
problem known as the Lur’e problem [20], [23], [25], [28]. 
The Popov and Circle criteria are sufficient frequency- 
domain conditions for absolute stability for the feedback 
interconnection of a continuous linear time-invariant system 
with a sector-bounded nonlinearity [6], [8], [15], [18], [27], 
[36]. Its discrete-time counterparts are known as the Tsypkin 
criterion [17], [24] and the Jury-Lee criterion [11]. Consider 
the nonlinear discrete-time dynamical system: 

x{k + 1) = Ax{k) + Bp{k), 

q{k) = Cx{k) + Dp{k), p{k) = -(j){q{k)), 

where A G B G C G D G 

and the nonlinear operator </> S $, where $ is a set of static 
functions that satisfy ^(0) = 0 and have some specified 
input-output characteristics. This paper considers several 


classes of diagonal nonlinear operators $ that are defined in 
Appendix m in which the individual nonlinearities are sector- 
bounded, slope-restricted, and/or odd. 

Suppose that the state-space system representation 
(A, B, C, D) is a minimal realization of the transfer func¬ 
tion G{z) — C{zl — A)~"^B -f D such that the triplet 
(A, B, C) is controllable and observableQ It is well known 
that the frequency-domain sufficient stability conditions— 
the Circle, Popov and Tsypkin criteria—for the absolute 
stability of the system ([T]i can be simply represented by linear 
matrix inequalities (LMIs) written in terms of the matrices 
{A, B, C, D} via the Positive Real Lemma and the KYP 
lemma [1], [24], [35]. In this paper the S-procedure converts 
the constraints on the nonlinearities into LMI conditions, 
in the same way as in the direct derivation of LMI-based 
tests for absolute stability [1], [22], [25], [34]11 A modified 
Lur’e-Postnikov function is proposed that produces less 
conservative conditions for global asymptotic stability. The 
resulting LMI problems are solved for several numerical 
examples using off-the-shelf software [7], [29], [19]. 


11. Sufficient Conditions for Asymptotic 
Stability 


A. Modified Lur’e-Postnikov Stability of a Discrete-time 
Lur’e System 

Due to the conservativeness inherent to quadratic Lya¬ 
punov functions when applied to nonlinear systems, this 
section considers the use of a modified Lur’e-Postnikov 
function that includes integral terms whose sign definiteness 
are implied by the sector-bounded property of the nonlinear 
operators ((p G To derive sufficient conditions for the 

origin of ([U to be globally asymptotically stable, consider 
the Lyapunov function: 


<ik,i 


V{xk) = xlPxk+ 2'^Qii 4>i{(T)da 




l^ia - (j)i{a)] da, (2) 
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^Alternatively, the controllability and observability conditions can be 
removed because strict inequalities will be used in the stability criteria (see 
Section 2 in [26]). Instead, the condition KerA fl KerS = {0}, which 
implies that the origin is the one equilibria, is assumed to hold. 

^Related results have been derived based on integral quadratic constraints 
(IQCs) [4], [12], [13], [21], [22], [26] for systems with repeated monotonic 
or slope-restricted nonlinearities in continuous-time cases. 


where 



Xk 


>11 

P 12 

Pl3 

A 

Xk = 

Pk 

, p^ = p^ 

^12 

P 22 

P 23 


qk_ 


^13 

^23 

P33_ 


Qii ^ 0 ; Qii ^ O 5 — I 5 • ■ • : 

and the subscript k indicates a sampling instance. Both pk 
and Qk are functions of the state variable vector Xk, and the 
above Lyapunov function is radially unbounded and positive 
for all nonzero Xk G K”. 

For P 12 = 0,Pi3 = 0 ,P 22 = 0,P23 = 0,P33 = 0,(3 = 
diag{(3ii} = 0, the Lyapunov function (|2]l reduces to the 
Lur’e-Postnikov function used in the derivation of the Popov 
criterion. All other Lyapunov functions considered in the 
Lur’e problem literature are subsets of (O. Some papers in 
the Lur’e stability literature that appear to have additional 
terms actually just introduce the S-procedure in a different 
way. The dependence of the Lyapunov function (|2| on the 
nonlinearities and the states provide additional degrees of 
freedom for reducing conservatism during stability analysis. 
The application of the S-procedure provides a standard 
way to consider more restrictive classes of nonlinearities 
without having modify the Lyapunov function. The degree 
of conservatism will be compared for sub-problems with 
different subsets of the Lyapunov function (|2]i, to assess the 
reduction in conservatism posed by each term. 

B. Lagrange Relaxations 

Yakubovich showed that the positiveness of a quadratic 
function /o(a;) in a constraint set expressed in terms of 
quadratic functions, say fi{x), i = can be im¬ 

plied by the relaxed form with (Lagrange) multipliers. The 
S-procedure is a special case of Lagrange relaxation in 
which the constraints are represented in terms of quadratic 
functions, so that the multipliers can be combined into an 
LMI inequality. The application of the S-procedure appears 
not only in control theory but also in the more generic 
optimization literature. For convenience, below is the form 
of the S-procedure used in the proofs of this paper. 

Lemma 1: (S-procedure for Quadratic Inequalities) For 
the Hermitian matrices 0^, i = 0,..., m, consider the two 
sets: 

(51) C0oC <0, VC e $ = {C e F"|C*0K < o, Vi = 
1 ,... ,m}; 

(52) > 0, i = 1,... , m such that 0o — YlT=i A®* < 0- 
The feasibility of (S2) implies (SI). 

C. Discrete-time Lur’e Systems with Slope-restricted Non- 
linearities 

While conic-sector bounded nonlinearities, which bound 
the global slope of the nonlinear functions have been 
heavily studied in the literature, such a description allows 
the local slope of the nonlinear functions </>( ) to vary 
arbitrarily from one time instance to another. That is, such 
nonlinearities do not imposes any local slope restriction. 
Most nonlinearities in practice have a local slope restriction, 
in which case a less conservative analysis condition may 


be achieved if these constraints on the nonlinearities are 
included in the analysis via the S-procedure. A local slope 
restriction also provides an upper bound on the change of the 
integral term in the Lyapunov function, provided that 0( ) is 
continuous almost everywhere (a.e.). 

Theorem 1: Consider a system of the form (ID with the 
memoryless nonlinearity (p £ H that is contin¬ 

uous almost everywhere. A sufficient condition for global 
asymptotic stability is the existence of a positive semidefinite 
matrix P — with a positive definite submatrix Pn = 
and diagonal positive semidefinite matrices Q, Q, T, T, 
N £ such that 


Gn 

Gi2 

Gi3 



p<T 

^12 

G 22 

C 23 

< 0 

(3) 

'-"13 

'-"23 

G33_ 




where 


Gn = + P13C + G^pS + G^P33G)A - Pn 

- P13G - G^pS - G^P33G + A^C^QiCA 

- C^Q^C, 

Gi 2 = A^(Pii + P13G + G^pS + C^P33C)B - P12 

- P13P - C^Pis - C^P 33 D - C^T 

+ {CA - CfQ + A^C^QS,CB p {CA - CfQ 

- C^Q^D + {CA - CfpN 

Gi 3 = A^Pi2 + A'^PisD + A'^C^Pls + A^C^P33D 

- A^C'^f - {CA - CfQ + A^C^QCD 

- {CA - CfpN 

C22 = B^{Pii P P13G P C^P^^ P C^P33C)B — P22 

- P23D - D'^P^3 - D^PssD - 2 C^T - TD 

- D^T P Q{CB -D)P {CB - DfQ - 
P B^C'^Q^CB - D^Q£_D 

P {CB - DfQ P Q{CB -D)- p-^Q - 2 N 
P Np{CB -D)P {CB - DfpN 
C23 = P'^Pi2 + B'^PisD P B'^C^PJs P B'^C'^PssD 
— B^C^P -f QD -f Qp ^ P p 
PQDP B^C^Q^D P 2 N 

- {CB - DfpN P NpD 

C33 = P22 + P23D P D^P^^ P D'^PssD - 2 C^f - TD 

- D'^f -QD- D'^Q - Qp-^ - p-^Q P D^Q^D 
- 2 N- NpD - D'^pN 

Proof: The difference in the Lyapunov function (O between 
the k P \ and k sampling instances is 


NV{Xk) = Ck {AlPAa - E^PEaKk 

r<lk+i,i 


rig 


4>i{(T)da 


( 4 ) 


rQk+i,i 
’ Qk,i 


[iicr - 0i(cr)] da, 








where 



Xk 


'A BO' 


7 0 o' 

a = 

Pk 

Pk+ 1 _ 

A - 

5 -^a — 

0 0/ 
GA CB D 

,Ea = 

0/0 
G P 0 


Slope restrictions on the nonlinearities place an upper bound 
on the first integral: 


2 '^Qu 


rQk+ 1 ,-. 


(j){a)da 


2=1 


riq 

< 2 ^ Q,, 

2=1 

= Cjc/iCfc, 


r(0fc+l,2 4^k,i){sik-\-l,i 9A;,2)'1 

I ~ 2/17 ~ / 


where C/i is given in Appendix HU Similarly, an upper bound 
can be derived on the second integral: 


riq 

2 ^Q^^ 1 

Mk+I 

1 

b 

da 



i=l 

(lk,i 





riq 


rlk+i.i 


.. fik+i, 

i 

= -2E 

Qii i 

1 cj>{a)da 

+ 2 

Eq« / 

^icda 

2 = 1 

J 

gk.i 


i=i ''gk.i 


2 = 1 

Qii ^ 


fk^i 

+ fk,i {qk-i-l, 

i Qk^i) 


riq 

+ ^ Qii^i [qk+i,i — ql,i\ 

i=l 

= ClU 2 Ck, 

where U 2 is given in Appendix HIl 

Since the (negative) feedback-connected nonlinearity is 
monotonic with slope restriction in addition to being [0, 
sector-bounded, i.e., (j) € H it can be shown 

that the following inequalities are satisfied at each sampling 
instance k and all indices i = 1,..., n,: 

^k,i Qk,i\ — 0; 

(5) 

{4^k+l,i 4^k,i) {qk+l,i — 0- 

(6) 

The following notations based on © are useful when apply¬ 
ing the S-procedure: 

n, 

^ (j)k,i = Cfe ‘S'lCfci C^) 

i=l 

n, 

^ ' ‘^^i'Pk+l,i[^i 'Pk+l,i qk+l,i] — Ck ^'iCkt (8) 

i=l 

where Si and S 2 are given in AppendixHB A similar notation 
based on the inequality (|6|l is: 

Uq 

i=l 

= cls3<:k, (9) 


where S 3 is given in Appendix HU 


Applying the S-procedure, if the LMI G = A^PAa — 
E'^PEa + Ui + U 2 — Si — S 2 — S 3 < 0 is feasible then 
AV{xk) < 0 is satisfied for the specific class of feedback- 
connected nonlinearities cj) G <1)^°’^^ All of introduced 

matrix (decision) variables are of compatible dimensions. □ 


Remark 1: The upper bounds on the integral terms in 
dUl are derived by considering the a.e. continuity of the 
nonlinearity In those upper bounds the slope-restricted 
properties of the nonlinear feedback are exploited, which 
gives sharper bounds than the upper bounds in [9], [10], 
[14]. 


D. Discrete-time Lur’e Systems with Slope-restricted and 
Odd Monotonic Nonlinearities 

Theorem [T] exploits more information on the nonlinear 
operator than the original Lur’e problem. This section derives 
a less conservative sufficient stability condition for the more 
restrictive class of nonlinearities that are odd monotonic, by 
the introduction of additional quadratic constraints. 

Theorem 2: Consider a system of the form ([TJ with the 
memoryless nonlinearity 0 G n fl ^odd, which 

is continuous almost everywhere. A sufficient condition for 
global asymptotic stability is the existence of a positive 
semidefinite matrix P = with a positive definite subma¬ 
trix Pii = Pii and diagonal positive semidefinite matrices 
Q, Q, T, f, N, L, and L G such that 


Gii 

Gi2 

Gi3 



'-^12 

G22 

G23 

< 0 

(10) 

'-^13 

'-^23 

G33_ 




where 

Gii = Gii 

Gi2 = Gi2 - (GA - CfL - {CA + CfL 

Gi 3 = Gi3 - {CA - CfL - {CA - CfL 

G22 = G22 - {CB - DfL - L{CB - D) 

- {CB + Dfi - L{CB + D) 

G23 = G23 - {CB - DfL -LD- 

- {CB - D)'^L -LD + 

G 33 = G 33 - d'^l -ld-ld- d'^l 
P roof: The difference in the Lyapunov function (|2]l be¬ 
tween the k -\- 1 and k sampling instances is the same as 
AV{xk) in dUi with the same definitions of Ck, Aa, and Ea- 
Sector-bounded, slope-restricted, odd-monotonic feedback- 
connected nonlinearities satisfy (l5]l and (|6]l at each sampling 
time k and all indices i = 1,... ,nq for cj) G n<l>ir’^^ D 
^odd and satisfy 


{qk+l.i qk,i')\fk+l,i P 4^k,i\ fky^ fk-\-l,i ^ 0 ; ( 11 ) 

{qk+l.i qk,i')\fk+l,i fk,i~^fk+l,i 

Si 

^ ‘^Qk-\-l.ifk-\-l,i- 


( 12 ) 










The following notations are motivated by ([IT]i and (Hall: 

riq 

^ ^ ^k+l,i {Qk+l,i Qk,i){^k+l,i T 

z=l 

= cISia, ( 13 ) 


and 

Uq 

^2L, 

i=l 

= cJ55a, 


Qk,i(^^k+l,i 4^k,i') Qk+l,i{^k+l,i 
^k,i'^4^k+l,i 



(14) 


where S 4 and S 5 are given in Appendix HJ 

If the LMl G = AlPAa - PEa + U 1 + U 2 -S 1 - 
S 2 — S 3 — S 4 — S 5 < 0 is feasible then AV{xk) < 0 is 
satisfied for the class of feedback-connected nonlinearities 
4> G n n ^odd- All of matrix (decision) variables 
introduced are of compatible dimensions. □ 


Remark 2: The inequality constraints over the odd mono¬ 
tonic nonlinearities (fTTIi and (fTSli can be found in [23], [31], 
[32]. The inequalities (ITTl l and (IT^ are written in terms 
of quadratic functions of C,k, so that the S-procedure is 
applied to combine the constraints with the negative definite 
condition over AV{xk)- 

Remark 3: (Frequency-domain inequalities) The well- 
known KYP lemma replaces an LMI with an equivalent 
frequency-domain inequality (FDI) in terms of the system 
transfer function [33]. The FDI provides a graphical tool 
to determine feasibility of the LMI for systems with one 
nonlinearity, i.e., Up = Uq = 1, such as used in the 
graphical implementation of the Popov stability criterion. 
It is straightforward to apply the KYP lemma to transform 
each LMI Q and (ITOI l into an FDI (not shown for space 
considerations). 

Theorem [T] and |2] described LMI feasibility problems. 
That is, each criterion considers whether the set {X G 
§"|T'(2f) < 0,-Y > 0} is empty or non-empty. The 
following lemma shows the equivalence between two sets 
of matrix decision variables. 

Lemma 2: (Feasibility considerations) The set of symmet¬ 
ric matrices affine in X, {X G S"|T'(-Y) < 0,Ar > 0}, 
is nonempty if and only if the set of symmetric matrices, 
{X G S”|T'(A:) < 0,X > 0} is nonempty. 

Proof: The (only if) part is obvious. To prove the (if) part, 
the LMIs in the set can be rewritten as: 

N 

= + (15) 

where X = X^ e N = n{n+l)/2, and the E, 

are of compatible dimension. Suppose that X^ G {X G 
S”|T'(X) < 0,X > 0}. With the definition X^ = X° + SI, 
X^ > 0 for any value of positive scalar <5 > 0. Now show 
that there exists a scalar J > 0 such that X^ is in the set 
{X G §"|T'(A:) < 0,X > 0}. Consider 

N 

= (16) 

i=l 


where Si G {0,<5} for each i — . ,N and the number 

of nonzero St is n. Since the eigenvalues of are 

continuous for 5 0, these eigenvalues approach the 

eigenvalues of as 5 —0. In particular, for sufficiently 

small (5 the eigenvalues of have negative real part so 

that T'(X^) < 0. □ 


Remark 4: (Computation for strict and non-strict LMIs) 
From the results of Lemma |2] any LMI solver that is 
guaranteed to converge for an LMI feasibility problem with 
strict inequalities will converge for the above LMI feasibility 
problems. 


III. Numerical Examples for Stability Analysis 


Numerical examples are provided to illustrate the analysis 
results in Section [III All LMI computations were performed 
with MATLAB’s LMI Control Toolbox. The sector bounds 
and slope restrictions in each element fi of the 
nonlinearities were taken to be the same. The maximum 
upper bound on the ^ was computed, where /i was taken to be 
linear in as shown in each example. A large enough value 
for ^ makes the LMI in Theorem [T] (or |2| infeasible. A small 
enough value for ^ makes the LMI feasible provided that 
the system is nominally stable, that is, the system without 
nonlinearities is stable. These two values provide upper and 
lower bounds for the value of for which the LMI switches 
from being feasible to being feasible. The maximum ^ for 
which the LMI in each criterion is feasible were computed 
by the bisection method. 

As discussed in Section m the proposed stability criteria 
were derived for a Lyapunov function that has more degrees 
of freedom than past stability criteria and so have the poten¬ 
tial for being less conservative. Table 1 shows that the criteria 
in Theorems 1 and 2 are less conservative for six numerical 
examples. As expected, the stability margins are larger when 
more information is provided on the nonlinearities connected 
with the LTI system, and Theorems [T] and |2] provide more 
accurate estimates of the upper bound on ^ or /r. Although 
Theorems [T] and |2] apply to semiproper and strictly proper 
systems, the numerical examples only include strictly proper 
systems to allow comparison to the other criteria in the 
literature, which assume D = 0. 

Example 1 is from [24] (p = 20: 


G{z) 


-0.5z-f 0.1 

(z2 -z-f 0.89) (z-f 0.1) 



TABLE I 

The Maximal LFpper Bound on the Sector Bound 


Stability Criterion 

Ex 1 

Ex 2 

Ex 3 

Ex 4 

Ex 5 

Ex 6 

Circle (</. e 

Tsypkin (</. e 

Haddad et al. ((/> € D 

1.0273 

0.18358 

0.21792 

2.91387 

0.03660 

0.03716 

1.0273 

0.18358 

0.21792 

2.91387 

0.03660 

0.03716 

1.0273 

0.18358 

0.21792 

2.91387 

0.03660 

0.03716 

Kapila et al. {(j) E H 

Park et al. (0 E H 

1.0273 

0.18358 

0.21792 

2.91387 

0.03660 

0.03716 

1.7252 

0.18358 

0.21792 

2.91387 

0.03660 

0.03716 

Theorem 1 {<j> e $7^’ D 

2.4475 

0.73082 

0.30203 

43.40412 

19.18289 

0.04613 

Theorem 2 (0 G D n ^odd) 

2.5576 

0.73082 

0.83686 

43.40412 

19.18289 

0.18975 


Example 2 has 5 states and 2 nonlinearities (p = 0: 


A = 


Example 6 has a wide range of pole locations including two 


B = 


C = 





poles near 

1 (u 

= 0: 


0.2948 

0 

0 0 O' 







0 

0.4568 

0 0 0 

A 

= diag{0.5359,0 

0 

0 0.0226 0 0 

? 



— 

0.0034 

0 

0 

0 0.3801 0 



'1 

0 

0 

o' 

0 

0 

0 0 -0.3270 



0 

1 

0 

0 

■-1.1878 

0.2341 ■ 




0 

0 

1 

0 

-2.2023 

0.0215 




0 

0 

0 

1 

0.9863 

-1.0039 

5 

B 

= 

1 

0 

0 

0 

-0.5186 

-0.9471 




0 

1 

0 

0 

0.3274 

-0.3744 




0 

0 

1 

0 

-1.1859 

1.4725 

-1.2173 -1.1283 -0.2611 



0 

0 

0 

1 

-1.0559 

0.0557 

-0.0412 -1.3493 0.9535 



1 

0 

0 

0 


Example 3 has a dense A-matrix and 2 nonlinearities (/i 

0 : 




'0.0469 

-0.3992 

-0.0835' 



A = 

0.3902 

-0.5363 

-0.2744 

5 



0.4378 

-1.3576 

0.4651 




'-0.5673 

-0.2785' 




B = 

0.1155 

-0.0649 





-2.1849 

-0.5976 




C = 

0.3587 

-1.0802 

-0.680: 

1 ' 


-1.3833 

-1.0677 

1.1497 


Example 4 has two poles at the same 

location (fj, = 20: 


'0.4030 

0 

0 


-0.2494' 

A = 

0 

-0.1502 

0 

,B = 

0.2542 


0 

0 

-0.1502 


-0.2036 

C = 

0.9894 

0.6649 

0.4339] . 



Example 5 has three poles at the same location (p = 2^ 


'0.4783 

0 0 

0 

- 


A = 


B = 


0 0.7871 0 0 

0 0 0.7871 1 

0 0 0 0.7871 

-1.5174 
1.2181 
0.2496 
-0.5181 

C = [0.8457 -2.0885 1.2190 0.1683] . 


C = 


1 1 0 0 0 0 0 0 0 
0 0 1 1 1 0 0 0 0 
0 0 0 0 0 1 1 0 0 
0 0 0 0 0 0 0 1 1 


The stability margins for Example 5 indicate that Theorem 
1 can reduce the conservatism compared to literature results 
by more than 53,200%. The stability margins for Example 6 
indicate that, for odd nonlinearities. Theorem 2 can reduce 
the conservatism by more than 400% compared to Theorem 
1 and even more when compared to literature results. 

IV. Conclusions 

This paper considers the analysis of global asymptotic 
stability for discrete-time Lur’e systems with nonlinearities 
of various types. The proposed Lyapunov function (|2]i is a 
superset of all past reported Lyapunov functions, and the 
new sufficient conditions for globally asymptotic stability 
are less conservative. Conservatism is further reduced for 
nonlinear operators that satisfy additional constraints in input 
and output behavior, such as being odd or having local slope 
restrictions. Numerical examples indicate that the reduction 
in conservatism can be many orders of magnitude. The 
results can be generalized to the analysis of L 2 -gain, root- 
mean-square gain, and decay rate in the standard way [1]. 
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Appendix I 

Definitions for Some Classes of Nonlinear Operators 


— {(t> ■ I (t>i{o') - ct] < 0 , Vct e R and for z = 1,..., n,}; 

— {4> '■ I 0 < \/0- ^ ij g and for i = 1,..., n,}; 

(7 — (T 

^odd — {(t> '■ R”’ —>■ K"’’ I </'i(CT) = —(j)i{—a), Vcr e K and for i = 1 ,... ,nq}. 

Appendix II 

Matrices for the Application of the S-Procedure in Theorems [T|and[2] 


t/i^ 


t/2 = 


{CA-CfQ 

Q{CB -D) + {CB - DfQ - Q^-i 

* 


A^C'^Q^CA - C'^Q^C 






-{ca-cYq 

QD + Q/r ^ 

-QD - D'^Q - Qn-\ 
A^C^Q^CB + {CA - C)^Q - C^Q^D 
B'^C'^Q^CB - D'^Q^D 
+{CB - DfQ + Q{CB -D)- 




* 




* 


0 

C^T o' 


0 

0 

A^C'^f 


* 

2^-^T + TD + D^T 0 


* 

0 

B^C^f 


* 

* 0 


* 

* 

2^-1 f + fD + D'^f_ 


0 

-{CA-CY 

N 



{CA - CYN 

^3^ 

* 

- {cb-dYn - 

- N{CB - 

D) 


-2V/X-1 + {CB - DY 


A^C^Q^D 

fi-^Q + QD + B^C^Q^D 
-fi-^Q + D^Q^D 


{CA - C^L 
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